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We present relativistic coupled-cluster (CC) calculations of energies, magnetic-dipole hyperfine 
constants, and electric-dipole transition amplitudes for low-lying states of atomic boron. The triva- 
lent boron atom is computationally treated as a monovalent system. We explore performance of the 
CC method at various approximations. Our most complete treatment involves singles, doubles and 
the leading valence triples. The calculations are done using several approximations in the coupled- 
cluster (CC) method. The results are within 0.2-0.4% of the energy benchmarks. The hyperfine 
constants are reproduced with 1-2% accuracy. 
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I. INTRODUCTION 

Atomic parity violation [1] (APV) provides powerful 
constraints on new physics beyond the standard model 
of elementary particles 2, 3]. While the experiments are 
carried out at low energies, the derived constraints are 
both unique and complementary to those obtained from 
direct searches for new physics with high-energy parti- 
cle colliders. For example, the latest APV results 0, H| 
set new mass bounds on exotic new-physics particles, Z' 
bosons, which are ubiquitous in competing extensions of 
the standard model. These APV bounds improve upon 
the earlier results of the Tevatron collider and cut out 
the lower-energy part of the discovery reach of the Large 
Hadron Collider. 

Interpretation of APV experiments in terms of elec- 
troweak observables requires input from atomic-structure 
calculations. APV is the field where focus on precision of 
both atomic experiment and theory is of importance. As 
we look at the entire body of experimental APV results 
with various atoms, we find that the most accurate mea- 
surements were carried out with 133 Cs atoms [1, 0|- The 
second best APV measurement was done with 205 T1 Q. 
Ideally, the uncertainty of calculations should match the 
experimental error bars. For Cs, the 0.27% uncertainty of 
relativistic many-body calculations J4,j5[ is better than 
the experimental 0.35% error bar [y, 0|. For Tl, how- 
ever, the situation is reversed. The experimental accu- 
racy here is about 1 % la, w hile the theoretical errors are 
estimated to 2.5-3% [9l-Fll| . While for Cs it took a little 
over a decade for the theory to match the experimental 
accuracy, for Tl, even after almost two decades the state 
of the theory remains inadequate. 

The goal of this paper is to start exploring the fea- 
sibility of transferring high-precision many-body tech- 
niques developed for Cs 0,[1] to the Tl atom. There is a 
crucial distinction between the Cs and Tl atomic struc- 
tures. Cs has a single 6si/ 2 valence electron outside the 
closed-shell core, while Tl nominally has three valence 
electrons (6s?, 2 6f>i/2 ground state configuration). Since 
6pi/2 is the only active electron involved in the measured 
6P1/2 — 6p 3 / 2 PNC amplitude [U, it is natural to wonder 
if the 6s^ 2 shell could be considered as a part of the core, 



thereby enabling reuse of the Cs techniques. The relevant 
figure of merit is the energy gap between the core and 
the valence subspaces. At the Dirac-Hartree-Fock (DHF) 
level, this gap for Cs (6si/ 2 — 5p 3 / 2 ) is 156479 cm -1 while 
in Tl, when treated as a monovalent system the gap be- 
tween the 6s!/2 and the 6pi/2 is about 1.5 times smaller. 



Considering the smallness of this gap, the interactions 
between the outer core-shell and the valence electrons 
have to be treated in non-perturbative fashion. 

Before proceeding to computationally-expensive 81- 
electron thallium atom, in this paper we focus on a proto- 
type atom, boron. Boron belongs to the same IIIA group 
of the periodic table but has only 5 electrons. The ground 
state configuration of boron reads ls^ 2 2s^ 2 2p 1 / 2 , with 
three valence electrons. Here, however, we have treated 
boron as a monovalent atom by placing the two lower en- 
ergy valence electrons (2s^ 2 ) in the core. Such a treat- 
ment greatly simplifies the underlying equations and cal- 
culations. Treating the boron as a monovalent system is 
also justified by the fact that the 2p electron is often the 
excited electron in optical transitions. 

Notice that in the non-relativistic Coulomb approxi- 
mation, the energies of the 2s and 2p electrons would 
be identical, again strongly suggesting that the perturba- 
tive monovalent approach would fail right from the onset. 
The mean-field effects lift this degeneracy resulting in the 
energy gap of 131292 cm -1 in the DHF approximation. 

Another compelling feature of boron is that as a five- 
electron system it lies at the applicability border of high- 
accuracy variational methods [12|. In the case of few- 
body systems, e.g., He and Li, the computational ac- 
curacies have been substantial (see, e.g., Refs. [l3|-[l5j].) 
Only very recently Bubin and Adamowicz (l2j have ex- 
tended the reach of the full-scale variational methods 
to boron. While such accurate variational calculations 
were carried out for boron, the problem with variational 
methods is the impracticality of extending them to even 
larger atoms, such as the 81-electron Tl. That leaves 
many-body methods as the best hope for accurate com- 
putations for Tl. Due to the strongly-correlated nature 
of these systems and the desired high accuracy, one has 
to employ non-perturbative (all-order) methods, where 
certain classes of diagrams are summed to all orders of 
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many-body perturbation theory in the residual Coulomb 
interaction between electrons. 

In this paper we employ arguably the most popular 
all-order method: the coupled-cluster (CC) method; this 
method was at the heart of high-accuracy APV calcula- 
tions for Cs [3, [H] and other heavy atoms . We will 
use an ab initio relativistic formulation. Qualitatively, 
various CC methods are distinguished by the maximum 
number of simultaneously excited electrons from the ref- 
erence DHF Slater determinant. Here we include single 
and double excitations of core electrons and single, dou- 
ble, and triple excitations of the core and valence elec- 
trons. We refer to this scheme as the CCSDvT method. 

With respect to the previous CC-type calculations for 
boron, highly-accurate non-relativistic calculations were 
reported in Ref. [TtJ ■ These authors included excitations 
of all five electrons (i.e., CC with single, double, triple, 
quadruple and quintuple excitations) and also used the 
Slater- type geminals at the CCSD levels to correct for the 
incompleteness of their necessarily very limited compu- 
tational basis. They then added relativistic corrections 
to their calculations in an ad hoc manner. These authors 
computed the ionization potential of the ground state of 
boron and other atoms. Less complete (CCSD trunca- 
tion level), albeit fully relativistic calculations were car- 
ried out in Ref. [l8j]. The focus of that work was on 
the convergence of iteration method for solving linearized 
CCSD (LCCSD) equations. Here we will use the conver- 
gence techniques developed in that work, but will employ 
a more sophisticated CCSDvT framework. 

This paper is organized as follows: In Section II we 
discuss the CCSDvT method and justify the use of the 
convergence method already implemented for LCCSD in 
Ref. [l8| . In Section III we present calculated energy lev- 
els of many of the low-lying states of born, their electric- 
dipole transition amplitudes and hyperfine constants. We 
also compare our results with other computational meth- 
ods and experimental data. Finally, in Section V we draw 
the conclusions. 



II. METHOD 
A. Coupled-cluster method and approximations 



DHF basis, ignoring a common energy offset, reads: 
H = H + G = y^£fiV[a|ai] + ]^^g l3k iN[a\a]aia k ] ■ 

i ijkl 

(1) 

Here Hq is the one-electron lowest order Hamiltonian and 
G is the residual Coulomb interaction. Ei is the single- 
particle DHF energy and gijki is the two-body Coulomb 
matrix element. N[...] indicates that the operators are 
in normal form with respect to the quasi-vacuum core 
state, |0 C ). Operators a' and a are respectively creation 
and annihilation operators. 

In our implementation of the coupled-cluster method the 
exact atomic wave function of a monovalent atom with 
the valence electron in state v is written as: 

i*„> = ni*(°>>, (2) 

where \^i 0) ) = a£|0 c ) is the lowest-order DHF wave func- 
tion. The wave operator ft maps the DHF solution onto 
the exact wave function, The ansatz for the wave 

operator is: 

fi = JV[exp(6)], (3) 

where C is called the cluster operator. For a system of 
N electrons it is expanded as: 

JV 
i 

where i indicates the number of excitations of core and 
valence electrons. For example, the operator G\ may be 
split into two classes of core and valence excitations: 

G\ = S c + S v , (5) 

Sc ^ ^ PmaQ j rn® , ai 
ma 

I 'l l v 

The coefficients p ma and p mv above are referred to as 
cluster amplitudes and are to be found. Here and else- 
where in this paper, the indices a, 6, ... are reserved for 
core orbitals, m, n, ... are designated to virtual or excited 
orbitals, v,w, ... indicate valence states and i, j, k, and 
I are arbitrary orbitals. We will subsequently equate C% 
with D and C3 with T . Explicitly, 



Here we briefly describe the coupled-cluster method 
and approximations used in our computations. We 
use coupled-cluster formalism for systems with one va- 
lence electron outside the closed-shell core. The reader 
may find more detailed descriptions of what follows in 
Refs.dHlj]. 

In our treatment of the atomic Hamiltonian for one 
valence electron systems, we employ the frozen-core 
Dirac-Hartree-Fock (DHF) potential. In the second- 
quantization notation, the atomic Hamiltonian in the 
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Expanding the exp(C) will lead to various powers and 
products of different cluster operators. Approximate so- 
lutions to the wave function are then found by keep- 
ing only a certain number of terms. As an example, 
the linearized coupled-cluster single-double (LCCSD) 
method keeps only the S and D terms: 



|tf„> = (1 + S + D)\$ v ) = OlccsdI^ 



(10) 



Another approximation, coupled-cluster single-double 2 
(CCSD2), keeps only up to the second-order terms in S 
and D and discards higher orders: 

\*v) = (1 + S + D + ^S 2 + X -B 2 + §D)\& V ) = Occsd 2 |^) 

(11) 

We used a combination of the above approximations with 
the addition of certain contributions of the triple term, 
T. Below we will discuss the approximate method used 
in more detail. 

In order to find the cluster-amplitudes, p's, one uses 
a set of generalized Bloch equations, see Ref. [2l[ for 
equations specific to a monovalent system. There are two 
such sets of equations for the core and valence states. The 
first set involves only the core-related operators. This set 
can be written as: 

( £v - ffo)^ COr 1^ 0) ) = ({QG^}conn.) COrO |*l 0) ) • (12) 

Here Q = 1 — | is a projection operator and e v 

is the Hartree-Fock energy of the valance electron. The 
subscript "conn." indicates that the retained Brueckner- 
Goldstone diagrams have no disconnected parts except 
for the valence lines. The above core equation does not 
depend on the valence state. 

The other set of Bloch equations is formulated for va- 
lence amplitudes, 



(e v + SE V - H )n™ l \¥^) = ({QGfi}co„„.) val |*r> 



\val|.-r.(0)\ 



(13) 



Here 6E V = {^„\QQ\^v' ) ) is the correlation energy. It 
should be noted that the right-hand side of the Eq. (TIB"]) 
contains both core and valence amplitudes. 

In order to solve the Bloch equations (fT2"j) and flT3"l) . 
one could devise an iterative approach. Doing this we 
obtain the recursive relation: 

(c + (SE V ) - H ) f>(" +1 >|*( )) = ({QGCl^}^.)^), 

(14) 

with = 1. In the above equation (SE V ) means that 
the 5E V should only be kept in the valence equation, (Tl3l) . 
Equations (TT2l to ([14")) are general and one can substitute 
different states in them. 

The core equations are solved by using the CCSD ap- 
proximation, Eq. (|11[) . In the present work, we used the 
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nonlinear terms in Eq. (fTTj) in addition to some lead- 
ing terms of the valence triples of Eq. |T7J All single and 
double contributions of the CCSD method (including the 
CCSD2 terms) are spelled out in detail in Ref. [22j and 
the triple terms, including some we have discarded, are 
shown graphically in Ref. (23|. Therefore, here we only 
discuss such contributions qualitatively. 

Both core and valence Bloch equations can be further 
separated into equations for single, double, and triple 
cluster amplitudes. For example, the to polo gical struc- 
ture of valence singles equation becomes |23| : 



■[H , S v ] + 8E V S V = CCSD2 + S v [f v 



(15) 



where the notation S v [T v ] stands for the effect of valence 
triples (T v ) on the right-hand side of valence singles (S v ) 
equation. It should be noted that in the equation above 
and what follows SE V = <5-EccSD2 + SE V [T V ] is the corre- 
lation energy. The equation for valence doubles reads: 



■[H ,D V ] + SE V D V « CCSD2 + D v [f v ] 



(16) 



Finally the approximation used for valence triple equa- 
tion reads: 

- [H , f v ] + SE v f v « f v [D c ] + f v [D v ] + f v [f v ] . (17) 

In the present work we use two different approximations 
for the right- hand-side of the triple equations, Eq flTTl) . 
If we keep only the two first terms on the right-hand- 
side of Eq. (HZ}, we call the method CCSDvT2, and if 
we additionally include the T^[T K ] terms, the method is 
called CCSDvT3. The reason for such a distinction is 
that the inclusion of T v [T v ] terms on the right-hand-side 
of Eq. (|17p . greatly increases the computational demands. 

We further simplify Eqs. p3| - pT|) by taking advantage 
of the spherical symmetry of atoms and analytically sum- 
ming over the magnetic quantum numbers. This results 
in a reduced form of cluster amplitudes and Coulomb in- 
tegrals. Wc use a different notation for such reduced am- 
plitudes. For example, the reduced valence amplitudes 
of singles are denoted as p(mv), doubles as pk{mnva) 
and triples as pk^^mnrvab). Here k\, k 2 are integer 
coupling momentum numbers and h is a half integer cou- 
pling angular momentum. As an example, the relation- 
ship between the ordinary and the reduced triples maybe 
represented as [22| 



J mnrvab 



E 

k\ k-2 h 



V ) 
K t 




— >— 

+ 


— S — 



Pk!k 2 h (mnrvab) 



™ (18) 

where the diagram subsumes various 3j symbols 19]. In- 
terested readers can find complete discussions of angular 
reduction in Refs. [l9| a nd J24T] . for single and double 
term reductions, and Ref. 22] for triple term reductions. 
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The angularly reduced equations have numerous terms 
on their right-hand sides and are composed of the sum- 
mation of products of cluster amplitudes, p's with each 
other and with two-body Coulomb matrix elements (g's 
in Eq. (JTJ).) The right-hand side terms can be found in 



Refs. [H, HE Hfil an d we will not be reproducing them 
here. After the angular reduction, the three valence 
Bloch equations (fT5 l) - ([T7l) can be represented in a gen- 
eral form as: 



(e m - e v + 5E v )p(mv) = ^ A s ■ p(nv) + ^ B s ■ p k (mnva) + ^ C s ■ Pk lk2 h(mnrvab) + ^ D s , (19a) 

n nak nrabkk'h 

{smn - £va + 5E v )p k (mnva) =y^ A D - p(nv) + y^B D - p k (mnva) + ) Cp ■ p kl k 2 h{mnrvab) + Dp, (19b) 



(£r. 



£ va b + SE v )pk 1 k 2 h(mnrvab) = ^ B T 



nak 



pk{mnva) 



jirabkk' h 



C T ■ pk^hjmnrvab) + } j D T , 



nrabkk' h 



(19c) 



where the amplitudes with tilde signs (p) are antisym- 
metrized combinations of the reduced amplitudes. In the 
above equations ery/. = + £j + £k and the A, B, C, 
and D's are constants with their subscripts denoting to 
which set of equations they belong, for example S for sin- 
gles etc. It must be emphasized that the right-hand sides 
of the above equations are linear in terms of valence clus- 
ter amplitudes and the constants include already known 
core cluster amplitudes and Coulomb matrix elements, 
g's. On the left-hand side, however, the terms are not 
linear in terms of valence cluster amplitudes due to the 
presence of SE V terms, which also depend on the valence 
cluster amplitudes. Despite this fact, in each iteration 
the SE V 's are taken from the previous iteration and sub- 
stituted to find a new set of cluster amplitudes and SE^s. 
In this sense the above equations are overall treated as if 
they were linear at each iteration. 



B. Convergence 



Simply iterating the CC recursive equation, Eq. ([M)) . 
at times proves to be insufficient as they could lead to 
slow convergence, non-convergence, or even convergence 
to nonphysical solutions. In Ref. we discussed two 
convergence methods developed in quantum chemistry 
(27j and their effectiveness in application to the LCCSD 
method. Here we apply one of the methods, the direct 
inversion of iterative space (DIIS) to our CCSDvT2 and 
CCSDvT3 methods. 



and 



To see how the DIIS is applied, we rewrite Eq. (|19a[) in 
a more streamlined fashion: 



p(mv) 
Pk(mnvb) 
Pkik 2 h(mnrvab) 




D 




Then the combined three sets of Eq. (|19a[) can be written 
as: 



Dt = a + At, 



(20) 



where A is a matrix including all the coefficients in 
front of cluster amplitudes on the right-hand side of 
Eq. (|19al) . The above equation can be solved iteratively 
by re- writing it as 



dm+l) 



D _1 (a + A-t (m >). 



(21) 



The iterative equation above is initiated by letting = 
on the right-hand side and finding and so on. 

The DIIS method works in two steps. First, a few 
iterative solutions of Eq. (|2Tj) are found. Next, a lin- 
ear combination of the said t' 1 ) is used to find the best 
next solution to the equation. For example, after accu- 
mulating m iteratively found solutions, v- 1 ', t^ 2 \..., t^ m \ 
the next best approximation can be found as their linear 
combination, 



■.(m+l) _ 



8=1 



ait 



W _ 



(22) 



Here Oi is the weight assigned to £w. In the case of 
the DIIS method, the a coefficients are determined by 
applying an error minimization scheme and solving the 
resulting system of linear equations [l8| 

T T (A - D) T a + T T (A - D) T (A - D)T<r = . (23) 

The new answer thus found, f( m+1 ) is then fed back to the 
right-hand side of Eq. (|2~T|) and the two steps are repeated 
until some parameter, i.e. the 6E V , stops changing (up to 
a specified accuracy) between two consecutive iterations. 
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C. Matrix elements 

After finding the cluster amplitudes and correlation en- 
ergies, we can calculate matrix elements of a one-particle 
operator 



Z = 



(24) 



In the present paper we refer to the results obtained by 
such a procedure as "scaling" . 

In the next section we will discuss the results obtained 
with different approximations and present valence ener- 
gies, matrix elements, and hyperfine constants for the 
boron atom. We will also compare our results with pre- 
vious experimental and theoretical results. 



where Zy is the single-particle matrix element. Notice 
that in deriving the CC equations one uses the interme- 
diate normalization scheme (*&t>|"^ 
elements then have the form, 



(0)\ 



1. The matrix 



M wv — 



y/(* w \* w )(* v \* v ) 



(25) 



As discussed in Ref. [25j, this matrix element could be 
separated into two parts leading to the expression 



y/[l + {5N W ) 

conn 

}[l + (SN v ) 

C 



(26) 



Here Z Q = (0 C \& z£l\O c ) and the remaining contributions 
of Z wv — (O c \a w Cl'l,ZCl v al\O c ) are encapsulated into Z\. 
In a similar way 



N w = (V W \V W ) = A H 

a = (o c |ntfi|o c ) = i, 



with SN W containing the rest of the contributions to the 
normalization. Zq contributions vanish for nonscalar op- 
erators Z, so we can ignore them here. In Refs. [25|, [28| 
the contributions to matrix elements in the LCCSD ap- 
proximation are explicitly listed. Ref. [29| further dis- 
cusses all the leading nonlinear contributions and the 
contributions of connected triple excitations to matrix 
elements. 

When the CC exponent of Eq.(j3|) is expanded in 
Eq. (|25[) , an infinite number of terms are produced. The 
resulting series may be partially summed so that it sub- 
sumes an infinite number of terms. This procedure is 
called "dressing" and is explicitly explained in Ref. (29j . 
In short, it is built on the expansion of the products 
of the CC-amplitudes into a sum of n-body insertions. 
Two types of insertions are considered: particle (hole) 
line insertion and two-particle (two-hole) random-phase- 
approximation-like insertion. It must be noted that this 
procedure is specialized for the CCSD approximation in 
monovalent systems. 

Due to the approximate nature of the CCSDvT 
method, certain correlation effects are lost. To partially 
account for the missing contributions in the calculation 
of the matrix elements, we correct the wave functions us- 
ing a semi-empirical procedure suggested in Ref. [30] | . In 
this procedure the valence singles are re-scaled by the ra- 
tio of experimental and theoretical correlation energies. 



III. RESULTS AND DISCUSSION 

So far we recapitulated various CC approximations. 
We also described how we apply the DIIS converging 
method to the iterative solutions of the LCCSD, CCSD 
and CCSDvT equations. In this section we present ab 
initio numerical calculations for properties of several low- 
lying levels of the boron atom. 

Atomic boron has three valence electrons, with the 
ground state configuration ls 2 2s 2 2p 1 . In our calcula- 
tions, we start by assigning the two electron 2s 2 valence 
orbitals of the ground state to the core orbitals. There- 
fore, we approximate the boron atom as a monovalent 
system. 

We employ the dual-kinetic balance B-spline basis set, 



obtained by solving the frozen core DHF equations [31 



This basis set numerically approximates a complete basis 
for the single-particle atomic states. The basis set was 
generated in a cavity of radius 40 a.u. and contains 40 
orbitals per partial wave for energies above the Dirac sea. 

The CC core equations are solved in the CCSD approx- 
imation. The core amplitudes are computed with partial 
waves summed up to and including the angular momen- 
tum I m ax = 6. The single and double core excitation co- 
efficients are then fed into the valence equations. The va- 
lence equations are initially solved in the LCCSD approx- 
imation. The resulting LCCSD wave functions are then 
used to initiate the nonlinear CCSD method. The valence 
CCSD result, in turn, becomes the reference point for the 
CCSDvT computations. In calculating the valence wave 
function, basis functions with l max — 6 are still used for 
singles and doubles. However, due to the computational 
expense of the triple terms, we employed a more limited 
basis set with l max = 3 for the triples. Initial calcula- 
tions with l max = 4 for all triple terms showed very little 
change in energies at a great expense in time. Therefore, 
to keep consistency of the data, we only present results 
with l max = 3 for the triple terms. Also, to keep the com- 
putational cost lower, we employed 35 out of 40 positive 
energy basis functions for the single and double terms, 
while reduced this number to 25 for the triples. Again, 
calculations of triples with 35 basis functions showed very 
little change in the outcome. We further introduced basis 
extrapolation corrections, l m ax — ► oo, which is discussed 
in subsection A. 

In the present work we use two different approxima- 
tions for triple equations. If we only keep the T V [D C ] and 
T V [D V ] terms on the right-hand-side of the Eq. (JTTJ), we 
call the method CCSDvT2. On the other hand, if we keep 
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all the terms in the Eq. (IT7|) we call the method CCS- 
DvT3 with f v [f v ] being the difference with CCSDvT2. 
The inclusion of T„[T„] terms on the right- hand- side of 
Eq. ([TT]) . greatly increase the computational cost of triple 
calculations. 

A brief analysis of the largest valence amplitudes (p's,) 
will follow in the next section (Sec. IIII A|) . Results for 
removal energies are presented in Sec. IIII Bl and for dipole 
matrix elements and magnetic-dipole hyperfine constants 
in Sec. iHTCl 



A. Largest contributions by cluster amplitude 

Before we present the results of the numerical calcula- 
tions, we would like to investigate the relative importance 
of individual valence amplitudes (p's) for each level of 
approximation. We run the numerical code, extract the 
five largest reduced valence amplitudes, and analyze how 
they change with each CC approximation. The results 
are shown in Figs. (Q~||3]). Here, we only analyze the va- 
lence single, double and triple amplitudes for the ground 
state (2pi/ 2 ) of boron. 

The valence states with the largest amplitudes remain 
for the most part the same from one approximation to 
another. For example, the CCSD approximation will ren- 
der the same five largest /5fc(mra;a)'s as the CCSDvT2 
approximation and so on. As is seen in Figs.([T][2]), in 
the CCSD approximation, the largest reduced valence 
amplitudes are lowered in value as compared to their 
LCCSD counterparts. Including triples in the CCSDvT2 
and CCSDvT3 approximations brings up these values to 
the region between the LCCSD and the CCSD approxi- 
mations. 

Based on perturbation theory, the cluster amplitudes 
are in general proportional to the ratio of Coulomb inte- 
grals (or their products) to energy differences between the 
orbitals. Therefore, we expect the orbitals with energies 
closer to the valence state under study to have the largest 
contributions to the valence amplitudes. However, as is 
evident from the figures, some of the larger contributions 
belong to states with high principal quantum numbers. 
For example, in Fig. [TJ the /5(10piy 2 , 2pi / 2 ) amplitude is 
the next largest after the p(3pi/2, 2pi/ 2 ). In order to un- 
derstand the unusually large contributions from highly 
excited orbitals, we plot the large component, P(r), of 
the 2p 1 / 2 , 2p 3 / 2 and 10px/ 2 orbitals in Fig.©, where r 
is the distance from nucleus. Examination of our DHF 
energies indicates that the continuum states start from 
the 5pi/2 orbital. The cavity size and its boundary con- 
ditions have the effect of compressing the continuum's 
10pi/2 wave function wave function toward the nucleus. 
In effect the Wpi/2 mimics the behavior of 2pi/ 2 in lower 
r regions, which in turn results in large Coulomb inte- 
grals between these states. The large contributions of 
other orbitals with high principal quantum numbers to 
the amplitudes is the result of such a cavity squeeze on 
the continuum wave functions. If we were to choose a dif- 
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FIG. 1: The five largest reduced single excitation coefficients 
p(mv) are compared in different CC approximations. Here 
v is the valence orbital 2p 1 ^ 2 and m is the excited orbital 
indicated on the x-axis of the graph. 
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FIG. 2: The five largest anti-symmetrized reduced double 
amplitudes, pk(rnnva), of the ground state of boron are com- 
pared in different CC approximations. Here v is the valence 
orbital 2p 1 / 2 and a is the core orbital 2s. Excited orbitals 
m and n and the angular momentum value k of the reduced 
amplitudes are shown on the x-axis. 



ferent cavity radius, some other continuum states would 
see a similar effect. 



B. Energies 

The computed energies of the 2p 1 / 2 , 2p 3 / 2 , 3si/ 2 , 
ipi/2, 3p 3 / 2 and 4s ^2 valence states of neutral boron 
are compiled in Table HI We compare our results with the 
National Institute of Standards and Technologies (NIST) 
recommended values [32| . In Table U we give a break- 
down of the contribution to the valence energies at each 
level of the CC approximation. For example, the energy 



7 



0.15 
0.14 
0.13 
0.12 
0.11 
0.1 
0.09 
0.08 
0.07 
0.06 





CCSDvT2 — — 




CCSDvT3 — — 











4 



(mnr) k 



FIG. 3: The five largest reduced triple amplitudes, 
Pk 1 k 2 h{mnrvab) , for the ground state of boron are compared 
in different CC approximations. Here v is valence orbital 
2pi/2 and a and b are the core orbital 2s. Excited orbitals m, 
n, and r and the coupling angular momentum values ki, 
plus the half integer h of the reduced triple amplitudes are 
shown on the x-axis. 



3 
CO 




2P1/2 

3p 1/2 
10p 1/2 



15 20 
r(a.u.) 



25 



30 



35 



FIG. 4: Large component, P, of the orbitals 2pi/2, 3pi/2, 
and 10pi/2 for boron is shown as a function of the distance 
from nucleus. The 10pi/2 wave function is squeezed by the 
cavity wall (40 a.u. here.) This makes the value of the 
p(10p 1( /2, 1p\i-i) amplitude comparable with p(3pi/2, 2p!/ 2 ) 
(Fig-O- 



difference between the DHF and the LCCSD methods is 
shown in the +LCCSD column, while the energy differ- 
ence between the LCCSD and the CCSD approximations 
is written under the +CCSD column heading and so on. 
The extrapolation values are calculated using the results 
in the LCCSD approximation. We first run the LCCSD 
code, limiting the angular momenta from l ma x = l to 6. 
Next we extrapolate the correlation energy difference be- 
tween each successive pairs of angular momentum num- 
bers by using the model discussed in Ref. [33[ . The values 



labeled "Total" are equivalent to the CCSDvT3 values or 
the addition of the contributions from all of the previous 
columns. Most of the correlation energy is recovered in 
the LCCSD method. Generally, as the complexity of the 
method grows, the respective correlation energy contri- 
butions become smaller. However, as can be seen from 
Table HI the energy contributions of the CCSDvT2 and 
CCSDvT3 methods are comparable to each other. The 
percentage difference between our final CCSDvT3 results 
and the NfST values are between 0.2% and 0.4%. 

In Table HH the same results are shown again, but this 
time the valence energies are tabulated with respect to 
the 2pi / 2 ground state. The results in Table |TT] are com- 
pared with the NIST recommended values, as well as two 
other theoretical calculations. The Gaussian-correlated 
(GC) method of Ref. [l2j is a non-relativistic approach, 
therefore, it docs not resolve the fine-structure split- 
ting for the p states. Multiconfiguration Hartree-Fock 
(MCHF) results of Ref. [34;] are calculated relativistically 
by using the Breit-Pauli Hamiltonian. The MCHF result 
for the 4s!/ 2 state, however, was not included in Ref. [34j ]. 
therefore, an older MCHF result [35] by the same group 
was used for comparison. 

The most accurate coupled-cluster calculation for the 
ground state of boron to date has been carried out by 
Kloppcr et al. 17]. The starting point of their method 
is known in the literature as the CCSD(F12) method. 
In their computational approach, the single-particle ba- 
sis sets are substantially truncated and are in effect in- 
complete. The incompleteness of the basis set is com- 
pensated by accounting for additional double excitations 
into Slater- type geminals (which is indicated by the F12 
qualifier) at the CCSD level of approximation. At this 
level of approximation our two methods should be equiv- 
alent to each other, since we saturate our basis sets by 
carrying out extrapolations to higher partial waves. How- 
ever, comparison of our and Ref. [17] results at individual 
levels of approximations is meaningless, since the start- 
ing point of the CC method, the independent-particle 
approximation, in our respective approaches is different. 
To obtain their high accuracy results, Ref. [17j include 
CC up to 5-fold connected excitations. In addition, they 
add relativistic corrections to their non-relativistic cal- 
culations results in an ad hoc manner. Their final result 
for boron's 2p 1 / 2 state is -66934.4 cm^ 1 . Nevertheless 
our CSDvT result without extrapolations is about 100 
cm -1 and with extrapolation 40 cm -1 off the experimen- 
tal value, while the final results of Ref. [13] recovers the 
experimental value. This indicates the important role of 
omitted higher-order terms (quadrupole and quintuple 
excitations) . 



C. Electric-dipole amplitudes and hyperfine 
constants 

After computing the wave functions as described in 
the previous section, we proceed to calculate the electric- 
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TABLE I: Contributions of each level of approximation to the valence energies of different B valence orbitals are shown and 
compared with the NIST recommended values [3(|. All energies are in cm" 1 units. 



B 


2pi/2 


2f>3/2 


3Si/ 2 


3^1/2 


3P3/2 


4*1/2 


DHF 


-60546.22 


-60528.30 


-25137.94 


-17258.14 


-17256.30 


-11368.93 


+LCCSD 


-6538.07 


-6537.01 


-1966.28 


-1244 


-1243.46 


-155.37 


+CCSD 


800.41 


800.14 


378.87 


228.44 


228.34 


-284.29 


+CCSDvT2 


-273.05 


-272.39 


-129.89 


-37.24 


-40.14 


-62.92 


+CCSDvT3 


-245.35 


-244.97 


-143.8 


-55.3 


-55.5 


-65.38 


Extrapolation 


-84.3 


-84.3 


-1.16 


-13.46 


-13.44 


-12.53 


Total 


-66886.58 


-66866.83 


-26993.76 


-18379.7 


-18380.5 


-11949.42 


NIST [32] 


-66928.04 


-66912.75 


-26888.35 


-18316.17 


-18314.39 


-11971.81 



TABLE II: Comparison of atomic energy levels of B, computed in different coupled cluster approximations, with NIST bench- 
mark values [3(| and two theoretical results: computational (multiconfiguration Hartree-Fock (MCHF) [3J] and Gaussian- 
Correlated (GC) [l2lp. All energies are in cm -1 units. GC method's results are nonrelativistic and do not resolve the fine- 
structure splitting, therefore there is only one GC value per total angular momentum, j. 



State 


LCCSD 


CCSD 


CCSDvT2 


CCSDvT3 


MCHF [34J 


GC [12] 


NIST [32J 


2pi/2 























2 P3/2 


18.98 


18.71 


19.33 


19.75 


15.39 




15.29 


3Si/ 2 


40069.65 


39648.11 


39791.23 


39892.82 


40005.27 


40048.20 


40039.69 


3^1/2 


48652.99 


48081.02 


48316.79 


48506.88 


49011.74 


48619.04 


48611.87 


3f3/2 


48655.39 


48083.32 


48316.19 


48506.08 


49013.49 




48613.65 


4Sx/2 


55631.76 


54547.06 


54757.15 


54937.16 


4642.89" 


55017.55 


55010.23 


"Ref. [H 



dipole transition amplitudes and the magnetic-dipolc 
hyperfine-structure (HFS) constants. We tabulate the 
results for the 2p 1 / 2 , 2p 3 / 2 , 3sx/2, and 4si/ 2 states of 
atomic boron. Here we discuss our results and compare 
them with the MCHF values of Refs. [34] and jH. 

In Table Hill we present the reduced electric dipole ma- 
trix elements, (3s||£)||2pj) and (4s||l)||2pj), computed in 
the length form [37j ■ Here, D is the electric-dipole oper- 
ator and j =1/2, 3/2. In the columns under the coupled- 
cluster heading, we list our results in order of the in- 
creasing complexity of the employed CC approximations. 
The corrections to the transition amplitudes, using the 
dressing and the scaling procedures, are in the next two 
columns. Our final results for the electric-dipole transi- 
tion amplitudes are written in the total column. These 
values are the addition of the corrections to the results 
obtained in the CCSDvT3 approximation. The MCHF 
method's result of Refs. 34 1 and [35| are tabulated in 



the last column. As can be seen from Table IIII1 our final 
results and the MCHF results are in a reasonable agree- 
ment with each other (about 0.3%.) 

In Ref. 34j, the stated values are in the linestrength 
form, Si, which is the reduced transition amplitude 
squared. The older Ref. |35| . contains non-relativistic re- 



sults for the reduced transition amplitudes between the 
4s and the ground state, 2p. In the LS coupling scheme, 
the relativistic and non-relativistic reduced matrix ele- 
ments, between s and p states, are related as 



\(n's 1/2 \\D\\np 1/2 )\ 
\(n's 1 /2\\D\\np 3/2 )\ 



-\(n's\\D\\np)\, 
~\(n' s\\b\\np}\. 



(27a) 



(27b) 



Therefore, we used the above relation to convert the non- 
relativistic results of Ref. [35[ to relativistic ones. 

A useful test of self-consistency of our results is to 
check the ratio of the reduced transition amplitudes, 
\(n's\\D\\npi/ 2 )/(n's\\D\\np 3 / 2 )\. This ratio should be 
equivalent to the ratio of the right-hand-side of Eqs. (|27a|) 
and (|27E|) . or l/y/2 ~ 0.7071. The ratio of the first pair 
of amplitudes in our different approximations is about 
0.7068 and for the second pair is 0.7076. These differ- 
ences can be explained by relativistic corrections. 

In Table HVl we compile the results of calculations of 
magnetic-dipole hyperfine-structure (HFS) constants, A, 
for the three states 2p 1 / 2 , 2p 3 / 2 , and 3si/2 of boron. Here 
we list the HFS constants for each level of the CC approx- 
imation. The corrections made to each HFS constant by 
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TABLE III: The reduced electric dipole transition-matrix elements (a||-D||6) in the length form in atomic units. 









Coupled Cluster 




Corrections 


Total 


MCHF 


\a) 


\b) 


LCCSD 


CCSD 


CCSDvT2 


CCSDvT3 


A (dressing) 


A (scaling) 






3s 


2pi/ 2 

2P3/2 


1.2454 
1.7618 


1.1772 
1.6656 


1.1623 
1.6444 


1.1542 
1.6330 


-0.0097 
-0.0143 


-0.0077 
-0.0105 


1.1368 
1.6082 


1.1345" 
1.6047° 


4s 


2^1/2 
2??3/2 


0.3736 
0.4791 


0.4420 
0.6251 


0.4498 
0.6359 


0.4480 
0.6331 


-0.0027 
-0.0037 


-0.0067 
-0.0093 


0.4386 
0.6201 


0A399 b 
0.6222 6 



the dressing, scaling, and extrapolation procedures follow 
the CCSDvT3 approximation's result. Our final results, 
the MCHF method's A constants [38[ as well as experi- 
mental results [39[ are shown at the bottom of Table IIVI 
Our 2pi/2 hyperfine constant is off by 2% from the ex- 
perimental result, while for the 2p 3 / 2 state the difference 
is about 1%. There are no experimental literature values 
for the 3s wo state of boron, so we compare our results 
with Rcf. [38] 's MCHF value. The percentage difference 
with MCHF calculation for the 3sx/2 orbital's HFS con- 
stant is about 0.3%. 



TABLE IV: Changes in valence magnetic-dipole hyperfine- 
structure constants of 2pi/2, 2P3/2, an d 3si/2 states of boron 
(in MHz) is shown for different CC approximations. A com- 
parison with experimental values of Ref. 39])and MCHF 
method's results of Ref. [3|| is presented at the bottom 
bracket. 



A 


2pi/ 2 


2P3/2 


3si/ 2 


DHF 


317.1 


63.3 


146.9 


LCCSD 


354.3 


87.0 


263.2 


CCSD 


358.6 


78.0 


235.3 


CCSDvT2 


364.1 


75.0 


240.2 


CCSDvT3 


368.0 


74.8 


242.8 


A (Scaling) 


0.7 


0.2 


-3.2 


A (Dressing) 


4.7 


-2.3 


-4.1 


A (Extrapolation) 


-0.2 


0.3 


1.5 


Final Result 


373.3 


72.7 


235.6 


MCHF [38] 


366.1 


73.24 


234.83 


Experimental [39] 


366.0765 


73.3470 





IV. CONCLUSION 



To reiterate, we examined the application of various 
coupled-cluster approximations for atomic boron. We 
treated the trivalent boron atom as a monovalent sys- 
tem, taking into account that the 2p valence electron is 
often the excited electron in optical transitions. We tabu- 
lated the results for a few valence energies, electric-dipole 
transition amplitudes, and magnetic-dipole hyperfine- 
structure constants of boron in the previous section. Fur- 
thermore, we compared our results with other computa- 
tional and experimental benchmarks. The results for the 
energies were found to be within 0.2% to 0.4% of the the 
NIST recommended values. The results for the electric- 
dipole transition amplitudes had about 0.3% difference 
with the MCHF benchmarks, while our HFS constants 
differed with the experimental values by 1% to 2%. Con- 
sidering that 0.1% accuracies are typical for true mono- 
valent systems (alkali-metal atoms,) the attained 1% ac- 
curacies for boron indicates deficiencies in treating it as a 
monovalent system. Indeed the comparison of our results 
with the more accurate CC computations of Ref. [13] 
shows that the way forward may be employing higher 
rank coupled-cluster amplitudes, quadruples and quintu- 
ples. 
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